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Abstract 



K^ , It has been known that Rb can be enhanced in the two Higgs dou- 

H I blet model if tan (3 is large. We point out that a similar enhancement 

in T{Z -^ T~^T~) is large enough to place a constraint on such a pos- 
sibility. We obtain a 95% CL upper bound ARh/Rh < 0.73% in this 
model for the MS mass mf,{mz) = 3.0 GeV. The 1996 world average 
is ARb/Rb = 0.97% ± 0.51%. We used the mb{mz) to determine the 
bottom Yukawa coupling instead of mb{mb) unlike in previous analy- 
ses, and also an improved experimental test of the lepton universality 
in Z decay, which made our results qualitatively different. 
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The Standard Model (SM) has been tested at an impressive accuracy 
by recent colhder experiments such as LEP, SLC, and Tevatron. Currently 
no serious conflict between the data and the SM is reported. However, the 
Z ^ bb branching fraction has been higher than the SM value at a few 
standard deviation level. The value of Rb = T{Z -^ hb)/T{Z —^ hadrons) re- 
ported in summer 1995 [1] was actually more than three standard deviations 
higher from the SM prediction and stimulated many theoretical and experi- 
mental efforts. The most recent measurements from SLC (0.2149 ± 0.0032) 
[2], ALEPH (0.2158 ± 0.0009) [3], DELPHI (0.2176 ± 0.0028 ± 0.0027) [4], 
and OPAL (0.2175 ± 0.0014 ± 0.0017) [5], however, are consistent with the 
SM {Rl^ = 0.2157 for mt = 175 GeV), and as a result, the 1996 world 
average value combining all the old and new data has come closer to the SM 
prediction {Rf, = 0.2178 ± 0.0011) [6]. Since older and newer measurements 
have different systematics, it is not clear at this stage whether it is appropri- 
ate to discard older measurements from the world average. In fact, DELPHI 
suggested to combine their older and newer measurements, while ALEPH did 
not. The central value may further evolve as newer methods are developed 
and the experimental inputs are updated. In view of this situation, it is 
premature to judge what the final outcome would be. Therefore, it is useful 
to investigate the consequence of various models on Rb in the light of other 
experimental constraints. 

It has been argued that the 1996 Rb value is actually more favored by 
various new physics scenarios than the 1995 one [7, 8, 9]. The 1995 average 
was hard to be explained by one-loop corrections to the Zbb vertex due to 
new particles. The 1996 average, on the other hand, is within the variation of 
Rb values predicted in the many new physics scenarios. This arises a renewed 
interest to check the consistency of various scenarios of high Rb with other 
existing experimental constraints. 

Currently there are two popular models which may lead to Rb values 
higher than that in the SM. One is the Minimal Supersymmetric Standard 
Model (MSSM) in the small tan/5 region, where the loop of a scalar top 
(mostly right-handed one) and a chargino (mostly higgsino-like one) can en- 
hance Rb which attracted many discussions [10]. The other is the loop of 
pseudo-scalar and scalar Higgs bosons in two Higgs doublet model (2HDM) 
which can also enhance Rb if tan /9 is large. The first scenario is now strongly 
constrained by recent limits on the chargino mass from LEP-II and the scalar 
top mass from DO. Still, one can accommodate a correction to Rb as large 



as ARb < 0.0017 [11]. It was argued that the latter scenario is in a conflict 
with the lack of an enhancement in the four-6 flnal states from Z decay [12] 
(see also [13]). The current experimental limit on four-6 or r^r^qq flnal state 
[14] is, however, not as stringent as estimated in [12], and there still remains 
a possibility that this scenario may enhance the Rb at a desirable level. It 
is the purpose of this letter to investigate whether this scenario is consistent 
with other existing constraints. 

The 2HDM has several motivations. First of all, it is the simplest exten- 
sion of the minimal standard model which has to be confronted by exper- 
iments. Second, the electroweak baryogenesis requires an extension of the 
minimal standard model to incorporate a CP-violation in the Higgs sector. 
The 2IIDM can naturally have CP-violating phases in its potential, and it is 
argued that it can create the value of cosmic baryon asymmetry as required 
by nucleosynthesis (for a review, see [15]). Third, the 2HDM may be a part 
of the MSSM. A general 2IIDM has a potential problem of flavor-changing 
neutral currents, which can be naturally avoided by either of the following 
two ways. The Type-I 2HDM lets only one of the Higgs doublets couple to 
quarks and leptons and hence the coupling matrix of the Higgs bosons can 
be simultaneously diagonalized as the mass matrix of quarks. There is no 
flavor- changing vertex of the Higgs bosons. In this case, however, a large 
tan/9 does not enhance the Higgs coupling to the 6-quark, and hence there 
is no interesting large contribution to Rb- On the other hand, the Type-II 
2HDM lets one of the Higgs bosons couple to the up-quarks while the other to 
the down-quarks and a large tan/9 can enhance the coupling to the 6-quark. 
In principle, either one of them can couple to the leptons. It is probably nat- 
ural, however, to assume that the one which couples to the down-quarks also 
couples to the leptons because they share the same weak isospin I3 = —1/2. 
It is indeed the case, for instance, in the MSSM. Then it is a natural question 
to ask whether a large tan/3 affects the phenomenology of the lepton sector. 

We find that there is a strong correlation between r{Z -^ t^t^)/T{Z -^ 
e^e~) and Rb in the Type-II 2HDM. From the observed lepton universality 
in Z decay, we find that Rb cannot deviate from the SM prediction by more 
than 0.73% at 95% confidence level almost independently from tan/9. This 
upper bound is in a confiict with the current central value of Rb- 

The consequence of the 2HDM on the Zhh and Zt^t~ vertices was studied 
already some time ago [16]. It was concluded that the correction to Zt^t~ 
vertex was too small to be observed compared to its size to Zhh vertex. We 



point out that two important changes should be made to this conclusion, 
however. The first is the improved accuracy in the experimental test of 
the lepton universality in Z decay, and the second is the running effect of 
the 6-quark Yukawa coupling between m\, and vriz which was not taken into 
account. A combination of a high accuracy in Y^ (defined below) and the 
running effect of my, can make the Zt^t~ vertex much more sensitive than 
previously thought. 

The lepton universality in Z decay is now tested at an extremely high 
accuracy. For later discussion, we parameterize the possible violation of the 
lepton universality by the following double ratio, 

_ r(Z ^ r+r-)/r(Z ^ r+r-)sM ,. 

^ v{z^in-)/v{z^in-)sM ' ^ ' 

Here, / refers to either e or /x assuming the universality among them, and 
the subscript SM to the standard model values. The advantage of us- 
ing Yt- is that many uncertainties cancel in the double ratio. Experimen- 
tally, the uncertainties in the luminosity measurement and overall width 
measurement nearly cancel between r and /. Theoretically, the top quark 
and Higgs boson masses enter the predictions of partial widths through 
oblique corrections, but they cancel in the ratio. The ratio of the SM val- 
ues V{Z -^ r+r^)sM/r(^ -^ l^l^)sM = 0.9977 is determined mostly by the 
kinematic factor /9^ = (1 — Am^/m^^Y^'^ to a very good approximation be- 
cause of the axial-coupling dominance in lepton couplings to the Z boson. 
Therefore, there is little theoretical ambiguity in the predicted ratio in the 
SM. We hence find that Y^- is the most useful quantity for our purpose. 

We derived the experimental value of Y^- from the Z line shape and lepton 
forward-backward asymmetries reported by the LEP Electroweak Working 
Group [6]. They quote measured values of m^, Tz, o"°. Re, Rfj,, Rt, ^fb' ^fb) 
Apg and their errors and correlation including new data from the 1995 energy 
scan. Here, Ri = T{Z — > hadrons)/r(Z -^ l^l~). Assuming the lepton 
universality for the first two generations, Ri = R^ = R^, Ap^ = Ap^ = Apg, 
and using x^ fit to the seven remaining parameters, we determine the ratio 

Ri T(Z^T+T-) 
The correlations among the seven parameters are fully taken into account. 



Normalizing it by the ratio in the standard model, we obtain 

K = 1.0008 ±0.0030. (3) 

The 95% confidence level upper bound-i- is Yr < 1.0057. Note that Ref. [6] 
quotes a somewhat weaker constraint on the ratio of electron and r couplings 
gvr/gve = 0.959 ± 0.046 and QAr/gAe = 1-0000 ± 0.0019, which correspond 
to 0.4% error in Yt-. However, the ratios do not assume e-/i universality and 
they tried to separate axial and vector couplings which are not necessary for 
our purpose. 

Second, we include the effect of the running of 6-quark Yukawa coupling 
between m^ and mz scales which was not taken into account in the previous 
analysis [16]. We take the following procedure. We first take the value 
of rrih = mh{mb) = 4.1-4.5 GeV in the MS scheme as summarized by the 
Particle Data Group [17]. To run the MS mass from rub to mz, we employ 
the renormalization group equation at two-loop level [18]. The numerical 
values oimb{mz) are shown in Fig. 1. The shaded region is the world average 
tts(^z) = 0.118 ± 0.003 [17].^ In this letter we take center values of m^ and 
as{mz), which correspond to mb(mz) = 3.0 GeV. We will discuss later how 
the results change for different choices of m^ and agimz). 

We do not go into the discussion of the MSSM Higgs sector in our letter. 
This is partly to simplify the analysis without many additional parameters in 
the MSSM. It is however mainly because a large contribution to Rb requires a 
large tan f3 with a light Higgs multiplets, which in turn implies a light charged 
Higgs boson in the context of the MSSM. Such a light charged Higgs boson is 
already strongly constrained by the CLEO measurement of the b —* s'-f rate. 
One needs to rely on a cancelation between the charged Higgs diagram (which 
always adds up with the standard model contribution) and the chargino loop 
which is also uncomfortably large. A light charged Higgs boson is further 
constrained also by B —>■ rvX [19] and r -^ /xz/z/ [20]. Even though there are 
viable regions in the parameter space [10], we do not pursue this direction 
further in this letter.^ 



"''The confidence level of the upper bound is determined by the one-sided Gaussian 
distribution. 

§Note that this average does not include the total Z hadronic width which is consistent 
with our spirit to allow Ri, to deviate from the SM prediction. 

^There is a possibility that the scalar tau loop may partially cancel the enhancement 
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Figure 1: Dependence of 6-quark mass at mz {rribijnz)) on asijnz). Two 
solid lines correspond to mfe(r77,fe)=4.1 GeV and 4.5 GeV, respectively. The 
shaded region represents the PDG average as{mz) = 0.118 ± 0.003. 

On the other hand, a general Type-II 2HDM allows a heavy charged Higgs 
boson naturally consistent with the 6 — > 57 constraint, while having neutral 
Higgs bosons as light as 50 GeV" . The phenomenological viability of such a 
light Higgs spectrum was recently also stressed in [23]. This is the parameter 
range of our interest in this letter. 

The 2HDM contains five physical Higgs bosons, two neutral CP even 
states h^ and if°, one CP odd state A^, and two charged states H^ and H~ . 

in T(Z —f T^T^) while keeping Rb large in a limited region of the parameter space. It is 
an interesting question whether such a cancelation is possible within the parameter space 
which sufficiently suppresses 6 — > 57 [21]. 

II The charged Higgs boson mass must be heavier than about 244 GeV due to the 
constraints from 6 — > S7 [22]. On the other hand, the electroweak p-parameter restricts 
the mass splitting among the Higgs bosons. For light h and A, we estimate the upper 
bound on the charged Higgs boson mass to be mH± ~ 200 GeV at 95% CL using the 
precision electroweak data given in [6] . Therefore, the 2HDM with light h and A is either 
only marginally consistent with these constraints, or requires some new physics to be 
consistent with the electroweak precision data. We simply assume in this letter that such 
new physics enters only oblique corrections and does not modify Zb or Zt^t^ vertices. 
We thank S.Kanemura for discussions on this point. 



We decompose two Higgs doublets as 

The mass eigenstates of the Higgs bosons are related to the weak eigenstates 
as follows: 

/i° = —rfi sin a + ri2 cos a, 
H^ = r]iCosa + r]2sma, 

for lighter (/i°) and heavier (H^) neutral CP-even Higgs bosons, 

A° = ^isin/5 + ^2Cos/5, 

for the neutral CP-odd Higgs boson, and 

H~ = H- sin/5 + {H+Y cos/5, 

for the charged Higgs boson. At this point, we have the following free pa- 
rameters: m/jo, m//o, 772^0, mH±, «, (3. Below, we take the limit of the heavy 
charged Higgs boson which is possible in the general 2HDM parameter space. 
We do so because the loop of the charged Higgs boson and the top quark 
reduces Rb and hence the positive contribution to Rb is maximized in this 
limit. Note also that the constraint from 6 — > S7 is naturally avoided in this 
limit as well. 

We evaluated the diagrams shown in Fig. 2, and found an agreement with 
previous calculations [16].** All the loop factors are common for b and r. The 
only differences are in the size of Yukawa coupling A/ and the couplings to 
the Z boson, gj^ = —\ — Qf sin^ 9w and g^^ = —Qf sin^ 6w for f = b,T. By 
parameterizing the loop factors from the first (second) diagram as e^A {^ff)-i 
the widths are corrected as 



AF; _ A^ 



> {gi? + {g^ 



2 

(5) 



R) 



^{gk- gL)^^hA - 4:gigieff 



**Wc ignored 6-quark and r-lepton masses in loop diagrams while we keep their Yukawa 
coupling constants finite. The vertex diagrams Z — > virtual {Zh'^ or H'^) — > bb are pro- 
portional to TObtan/3 but lack the tan^ /3 enhancement. For large values of tan/3 we are 
interested in, this contribution is subdominant and the omission of mh in the diagram is 
justified. The diagrams with Nambu-Goldstone bosons in the R^ gauge are the same as 
in the SM and hence are not enhanced for large tan/3. 
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Figure 2: Feynman diagrams with the neutral Higgs bosons in the 2HDM 
which contribute to the partial widths Z -^ bb and t^t^ and have the tan^ (3 
enhancement. 



Both thA and e// turn out to be positive. For both b and r, gj^ — gi = 1/2 and 
hence the first term in the bracket is simply etA/'^- On the other hand, the 
second term is 0.25e// for r and 0.13e// for 6, and hence is less important for 
b. In order to simplify the discussion we take a = (3 where the H^ is almost 
decoupled from leptons and down-quarks. In this case h^ has the enhanced 
Yukawa coupling for large tan (3 while H^ is the SM Higgs boson. The value 
of iriHO is irrelevant to the following discussion. ^^ One can easily generalize 
the analysis for a ^ (3, but the strong correlation remains the same and 
hence the final conclusion as well. The discussion then depends only on the 
following three parameters: tan/9, nij^o, and m/jO. 

The correction to Y^ from the neutral Higgs loops is shown in Fig. 3. 
This is a contour plot of the AYr induced from the neutral Higgs loops 
in the {771^0 ^niAo) plane with tan/? = 70. The shaded region is excluded 
by the negative direct search for Z -^ A^ h^ {rriho + tuao > mz)- The 
present experimental bound AY^ < 0.0057 is shown by the solid line, and is 
competitive with the constraint from the direct search. The excluded region 
from Yr is wider for larger tan/? because the correction is proportional to 
tan^ (3. 

The correction to Rb from the neutral Higgs loops is shown in Fig. 4, 
as a contour plot of the ARf,/Rh induced from the neutral Higgs loops with 

^^^^Note that there is a dependence on the mass of the standard model Higgs boson at 
0.2% level when one predicts partial widths theoretically. However this dependence comes 
through oblique corrections and cancels in the ratios Rh and Yr. 
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Figure 3: A contour plot of the AYr/Yr induced by the neutral Higgs bosons 
on the (m/jO,m^o) plane with tan/? = 70. The shaded region is excluded by 
the negative search for Z ^^ A^ h^ {f^h" + f^A" > f^z)- 
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Figure 4: A contour plot of the ARh/Rh induced from the neutral Higgs loops 
on the {mho,mAo) plane with tan/5 = 70. The shaded region is excluded by 
the negative search for Z ^ A^ h^ {^^ho + itia" > iti^z)- The solid line 
represents the bound from AY^- < 0.57%. 



tan/3 = 70. The solid line represents the 95% CL upper bound AYr < 0.57%. 
Recall that the central value of the 1996 world average is ARb/Rb = 0.97%. 

As can be seen from Figs. 3 and 4, both ARb/Rb and AY^- become maxi- 
mum in a region with tti^o = 171^0. This is particularly true for Rb because eff 
is less important and ehA is maximized when m^o = mf^o. In order to obtain 
the most conservative constraint, we take m^o = 771^0 which maximizes Rb 
while keeping Yj. small. In Fig. 5 we show the correlation between ARb/Rb 
and AYr with tan/9 = 50, 70, 90, assuming ttiao = iriho. Marks in each lines 
correspond to tjiao = rriho = 50, 60, 70, 80, 90, 100, 150 GeV. The shaded re- 
gion is the current experimental value of AYr at one standard deviation, and 
the solid line is the 95% upper bound AYr < 0.57%. From this correlation, 
the current bound on AYr constrains ARb/Rb < 0.73%. It is interesting that 
the 1996 world average ARb/Rb = 0.97% is well beyond the 95% CL upper 
bound. 

Finally, we would like to discuss dependence of the upper bound on Rb 
on the 6-quark mass and asimz)- The upper bound on Rb is proportional 
to m1{mz)-, and then larger as{mz) or smaller rrib gives a more stringent 
constraint. If we take q;s(t?2^)=0.121 (0.115) with mb=4.3GeV, it can be 
found from Fig. 1 that the upper bound on ARb/Rb becomes 0.71% (0.75%). 
A lattice calculation on T spectroscopy gives Wib = 4.0 d= O.lGeV [24], and a 
recent QCD sum rule analysis using more recent data of the electronic partial 
width of T(9460) in Ref. [25] favors a smaller value 4.1 GeV than that in 
Ref. [18]. If we use a relatively small value rTib =4.1GeV, the upper bound is 
reduced to be 0.65%. On the other hand, the Heavy Quark Effective Theory 
gives a lower bound on rfib {fUb > 4.26GeV in Ref. [26] and frib > 4.2GeV in 
Ref. [27]). If mfe=4.5GeV, the upper bound weakens to 0.81%. 

In summary, we pointed out that there is a tan^ f3 enhancement in r{Z -^ 
T^T^) induced by the neutral Higgs boson loops in the Type-II 2HDM, when- 
ever there is a similar enhancement to Rb- We found that the current experi- 
mental upper bound Yr < 1.0057 at the 95% confidence level places an upper 
bound ARb/Rb < 0.73% (for m^ = 4.3 GeV and asimz) = 0.118), which can 
be compared to the 1996 world average: ARb/Rb = 0.97% ± 0.51%. 
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